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In this lecture, we investigate the relationship between total dual integrality and integrality of 
polytopes. We then use a theorem on total dual integrality to provide a new proof of the Tutte-Berge 
formula. 

1 Total Dual Integrality 

Consider the linear program defined as 

max c T x 

s.t. Ax < b ^ > 

where A and b are rational and the associate dual program 

min y T b 

s.t. A T y = c (2) 

y>o 

Definition 1 The system of inequalities by Ax < b is Total Dual Integral or TD I if for all integral 
vectors c the dual program has an integral solution whenever the optimal value is finite. 

The main result for today is 

Theorem 1 If Ax < b is TDI and b is integral then P — {x : Ax < b} is integral 

Proof: We proceed by contradiction. Consider a vertex x* of P such that x* £ 7L. We can 
construct an integral c such that x* is the optimal solution corresponding to c by picking a rational 
c in the optimal cone of x* and scaling. Consider c = c+ |ej where q is an integer. Since the cone 
is full dimensional, c will still be in the optimality cone of x* for q sufficiently large. Now it follows 
that qc= qc + ej and thus (qc) T x* — (qc) T x* = x* £ Z. This means that either (qc) T x* or (qc) T x* 
are not integral which contradicts the assumption of total dual integrality. □ 
Note that the converse doesn't generally hold. We can have Ax < b integral with b an integral 
vector, but the system is not TDI. 

1.1 Total Unimodularity 

As an aside, we can consider an alternate condition which guarantees integrality. 

Definition 2 A matrix A is totally unimodular (TUM) if for any square submatrix A', det A' e 
{-1,0,1}. 

The following propositions hold for TUM matrices. 
Proposition 2 If A is totally unimodular then for all integral vectors b, Ax < b is integral. 

This differs from Total Dual Integrality where the integrality was dependent on both A and b. 
Proposition 3 If A is totally unimodular, Ax <b is total dual integral for any integral vector b. 
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For the case of non-bipartite matching, A is not totally unimodular. Indeed, for the three cycle 
with edges ei = (1, 2), e 2 = (2, 3), e 3 = (1,3), the matrix of x(S(v)) < 1 is given by 



A = 



1 1 
1 1 
1 1 



(3) 



and det A = 2. 

Often we can find subsystems of inequalities which define specific solutions of interest and are 
totally unimodular. Such a technique is employed in Lovasz's proof of the Lucchesi- Younger theorem 
on dicuts. 

1.2 An alternate proof of Theorem 1 

We begin by proving a theorem of Kroneckcr. 

Theorem 4 (Kronecker Approximation Theorem (1884) ) Ax = b has an integral solution 
if and only if y T b is an integer whenever y T A is an integral vector. 

Proof: To prove the forward implication, take an integral solution x*. Then y 1 Ax* = y T b and if 
y T A is integral then y T b must be an integer. 

To prove the converse, first note that there must be some solution to the system of equations; 
otherwise there would be a solution to y T A = with y T b ^ and by scaling y, we can get y T b £ Z. 
For the remainder, we will consider only a full row rank part of A. 

We proceed by introducing operations on the matrix A which preserve integrality. Let the jth 
column of A be denoted by aj. First note that exchanging two columns of A preserves both the 
existence of an integral solution of Ax = b and the property that y T b G Z whenever y T A E Z. Second, 
note that we can add any integral multiple of one column to another column and still preserve the 
assumptions. Indeed, for A e Z, if Ax = b, construct the matrix A' with columns identical to those 
of A but with a'i = ai + Xaj . Let x' be a vector with x' k — Xk except for x^ = Xj — Xxi . Then it is 
clear that A'x' = b and x' is integral (whenever x is). Conversely if A'x' = b, we can define x by 
x k = x' k except for xj = x 1 - + \x\ and x is integral and satisfies Ax = b. The preservation of the 
second assumption is proved similarly. 

Using these elementary operations, we can transform A into the form 



with B lower triangular as follows. For the first row, we can pair any two nonzero entries and 
compute their gcd using Euclid's algorithm 



since these operations are elementary, we can perform them on the columns and reduce the first row 
to one nonzero entry. We can then put this column as column 1 and proceed to the next row leaving 
column 1 fixed. Proceeding in this manner results in the desired form for A' . 

Now observe that B is nonsingular because we have assumed that A has full row rank. B~ x A' = 
[I ] and hence B~ x b must be integral (since every row of B^ 1 is a possible candidate for y T ). 
Since 



we have found an integral solution to the system A'x = b and this completes the proof. 



A' = [ B ] 



(4) 




(5) 



□ 
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Corollary 5 P = {x : Ax < b} is integral if and only if each supporting hyperplane contains and 
integral vector. 

Proof: The forward implication is immediate because every supporting hyperplane contains a 
vertex of P. For the converse, suppose x* is a non-integral vertex of P. x* is a unique solution 
of a subsystem Ax = b and by the Kronecker approximation theorem, there exists a vector y such 
that y T b is non-integral and y T A is integral. By adding an integral constant to the components 
of y, we can assume that y is nonnegative. Let c — A T y and a = y T b. Then c T x — a is a 
supporting hyperplane (teh fact that c T x < a is valic follows from the nonnegativity of y) and c T x 
is non-integral for all integral x which is a contradiction. □ 

This results in a new proof for Theorem 1. 
Proof of Theorem 1: If Ax < b is TDI and b is integral, pick an integral c such that a and Cj 
are relatively prime for i ^ j. By linear programming duality, maxc T x such that Ax < b will be an 
integer a and c T x — a will be a supporting hyperplane. 

Since the entries of c are relatively prime, we can find an integral vector x contained in the 
supporting hyperplane. (Indeed, it can be shown easily by induction on n that if the gcd of the 
entries of c is g then there is an integral solution to c T x — g.) Therefore, we conclude that Ax < b 
is integral. □ 

2 Back to matchings 

Given a graph G and a matching M define a vector x M £ ^-' E ' as 

M f 1 eeM 

e [0 otherwise v ' 

The Matching Polytope is the convex hull of all such incidence vectors. 
Consider the polytope P defined by the inequalities. 

x(8(v)) < 1 Vw e V 
\U\ 



x{E{U)) < 



_ 2 

x > 



VU € Todd (8) 



where V dd denotes the odd cardinality subsets. 

Edmonds proved in 1965 that P was indeed the matching polytope. Cunningham and Marsh in 
1978 proved Edmonds result by showing that P was TDI. Indeed this immediately implies that P 
is the matching polytope because all vertices of P would be integers, and any valid integer solution 
of P is a matching. 

Explicitly we have 

Theorem 6 (Cunningham- Marsh) For all w e l) E \ there exist integral vectors y and z such 
that the maximum weight of any matching is equal to 



mm^2y v + ^ 

vev sev od 



vev sev odd 

y > 0, z > 



5-3 



The proof of this theorem is in Schrijver's book. There are actually two proofs in the book, one 
that assumes the knowledge of the matching polytope, the other that's self-contained. 

We will now show that the Cunningham-Marsh theorem implies the Tutte-Berge formula in the 
cardinality case (w e = 1 for all e). 

Theorem 7 (Tutte-Berge) 



KG) = min - (\U\ + \V\ + o(G - U)) 
ucv 2 

Proof: Recall from lecture 1 that "<" was immediate. 

For any solution of the Cunningham-Marsh dual problem, it is clear that y v and zs are at most 
1 as w e = 1 for all edges. Furthermore, for all v e V cither y v = 1 or z$ = 1 for some odd set S 
containing v. 

Suppose zs and zt are such that SflT^fl. If S U T is an odd set, we can set z$ = zt = and 
zsut — 1 because 



\S\ 



+ 



\S\ + \T\ 



1 > 



\SUT\ + 1 



1 



ISUTI 



(9) 



and this assignment would reduce the cost function. If S U T is an even set, take j <E SUT and set 
z suT-{j} — 1 an d Vj = 1- This will also never increase the cost function. Therefore, we conclude 
that for an optimal solution, the sets {S G V dd ■ z$ = 1} are not overlapping. 

Let U = {v e V : y v = 1} and W = {S £ V od d ■ z s = 1}. If v e U and v e S with S eW, then 
we can remove v and an additional vertex u from S and let y u = 1 and this gives another feasible 
solution without increasing the cost function. Thus we can assume that U and all the sets S in W 
are disjoint. This implies that there cannot be any edges between the sets with zs = 1, which means 
that \W\ = o(G — U). Therefore we have shown 



vev sev od 



zs 



\u\ + vtim-\\w\ 

l -(\V\ + \U\-o{G-U)) 



(10) 



□ 
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